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SUMMARY 


The  method  Is  based  on  a  cornparlson  betvTeen  the  results  of  double 
integration  of  accelerometer  output  and  positions  as  observed  by  grovind  based 
instrumentation.  The  effects  of  errors  in  ancillary  equipment  are  discussed 
and  tables  of  accelerometer  errors  have  been  computed  for  a  typical  Black 
Knight  trajectory.  It  appears  that  the  accuracy  of  the  ground  eqTiipment  at 
present  available  at  Uoomera  is  barely  adequate  to  meet  the  requirements  of 
the  experiment. 
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1  INTRODUCTION 

The  difference  between  true  position  and  that  found  by  double  integra¬ 
tion  of  accelercneter  output  may  be  expressed  as  a  fionction  of  accelerometer 
misalignment  and  platform  drift.  An  over-determined  set  of  equations,  based 
on  this  relationship  between  integrated  values  of  recorded  acceleration  and 
position  observed  by  ground  instrumentation  for  various  times  during  flight, 
gives  a  means  of  evaluating  errors  in  the  accelerometers  themselves.  Further, 
if  the  r.m.  s.  errors  of  observation  are  known  it  is  possible  to  predict  the 
accuracy  with  which  the  accelerometer  errors  may  be  determined. 

This  has  been  done,  using  a  basic  Black  Knight  acceleration  curve.  An 
estimate  of  the  accuracy  with  which  the  accelerometer  bias  and  scale  errors 
can  be  calculated  is  made  and  tables  have  been  dravm  up  showing  the  weights 
which  should  be  given  to  the  probable  experimental  errors,  both  randan  and 
systematic,  in  the  ancillary  equipment. 

2  CCgJPONENTS  OF  ERRCE  IN  THB  ACOF.IFROMETERS 

Errors  which  are  anticipated  in  the  accelerometers  are  given  below  with 
the  notation  to  be  used  subsequently. 

(1)  Bias  k^. 

(2)  Constant  scale  error  k^. 

(3)  Scale  error  varying  v/ith  acceleration  k^x,  v/here  is  constant. 

(4)  Cross  coupling  errors. 

(5)  Misalignment  errors. 

For  the  Kearfott  accelerometers  with  v^hich  the  equipment  is  concerned 
initially,  the  expected  values  are 

k^  la  2  X  10  ^g  6*4  X  10  ^  ft  sec  ^ 

k^  2  X  10  ^ 

k2  ^  2  X  10“^  ft”*”  sec^  . 

The  cross  coupling  error  when  measured  in  laboratory  tests  was  small 
compared  with  these  and  has  been  ignored.  In  theory  the  axes  of  the  accelero¬ 
meters  will  coincide  with  the  fixed  space  axes  defined  at  launch.  In  practice 
there  will  be  errors  in  lining  them  up  but  the  discrepancies  will  be  measured 
and  should  be  known  to  5  seconds  of  arc. 

3  NOTATION 


The  axes  used  are  fixed  space  axes  and  defined  as  0,XIZ  whore  0  is  the 
centre  of  the  table  and  OZ  is  the  local  vertical  at  launch,  OY  the  direction 
of  the  line  of  fire  and  OX  completes  the  right  handed  set  of  orthogonal  axes* 
It  v^ill  be  assimaed  that  the  accelerometers  are  at  the  centre  of  gravity  of  the 
missile;  also  that  all  positions  determined  fran  ground  instrumentation  will 
have  been  corrected  for  the  rotation  and  ciurvature  of  the  earth  and  refer  to 
the  fixed  space  axes. 
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X 

X 

acc. 

•  • 

Sxt 

®yt 

Szt 

''gt 

"gt 


A 

2  V  r 


true  acceleration  in  the  direction  OX  at  time  t 
acceleration  experienced  by  the  x  accelerometer  at  time  t 
recorded  accelerometer  output  from  the  x  accelerometer  at  time  t 

gravity  components  at  time  t 


position  fran  groiand  instrumentation  at  time  t 


angles  betvfeen  the  axis  of  the  x  accelerometer  and  OX,  OY,  OZ 
at  launch  v/here  0  is  small 


4t 


t  «m  ,n 
x’  x'  X 

,m  ,n  I 

y  y  y  ^ 

■S  ,m  ,n 
z'  z’  z 


the  components  of  table  drift  about  the  fixed  axes  QPCYZ  at 
time  t.  (j>  small 


the  direction  cosines  of  the  accelerometer  axes  at  launch  so 
that 


m  :C:  6 

X  xy 


n  0 
X  xz 


^x 

^1x 

k2x^ 

^xt  ^ 


bias  error  in  the  x  accelerometer 
constant  scale  error  in  the  x  accelerometer 

scale  error  varying  with  acceleration,  kg  constant 

+  k .  X  +  X  ^ 

Ox  1x  acc.  2x  acc. 

X 

acc. 


T 


A 


ran 


in 


a 


m 


a  ♦ 

m 


a,  0 


position  and  velocity  error  components  along  the  x  axis  at  time  T 
a  sampling  interval 

a  rectangvilar  matrix  of  dimensions  m  x  n 
a  diagonal  matrix  of  dimensions  m 
a  column  vector  of  dimension  m 
a  row  vector  of  dimension  m 

the  angular  displacements  of  a  star  in  the  field  of  view  of  a 
television  camera 


m  ,n 


the  direction  cosines  of  the  line  of  sight  of  a  star  relative  to 
the  camera  axes 
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f  the  focal  length  of  the  camera 

e  a  random  error 

T+t'  T+t'  u 

j  xdt^  s  f  du  f  x(v)  dv 
T‘  T*^  T 

!».  METHOD  OF  EVALUATING  THE  COMPONENTS  OF  ERROR  IN  THE  ACCELEROMETERS 


The  method  described  for  determining  errors  in  the  aocelerometers  has 
been  developed  from  that  suggested  by  Chester^  which  was  found  to  be 
impracticable  when  applied  to  a  standard  Black  Knight  trajectory  during 
thrust.  All  accelerometer  errors,  drift  angles  and  misalignment  angles  are 
supposed  small  and  second  order  terms  are  neglected. 

At  time  t  the  direction  cosines  of  the  X  accelerometer  axis  with 
respect  to  0,XYZ  will  be 


6  .■&  +  m  -  6  .n 

^zt  X  X  ^xt  X 


~‘<f>  +  n  ) 

’^yt  X  ’^xt  X  X 


or  to  the  first  order. 


(1 ,  ^  ,  +  0  , 

'  '  ^zt  xy 
accelerometer  will  be 


+  ®xz^  acceleration  expeid.enced  by  the 


Now 


X  =  x  +  g  +  (^.+0  )(y  +  g)+(-^6_^.  +  0  )(z  +  g) 

aco.  ‘=’x  ^’^zt  xy'  ^  ^yt  xz'  ^  ^z' 


X  =  X  +k  +k  k  +k„X 
r  acc.  Ox  1x  acc.  2x  aco, 


and  this  can  be  expressed  as 


X  =  X  +  k  .  X 

r  acc.  xt  acc. 


where 


Ox  '  ■""Ix'^aoc,  '^x  acc. 


k^  +  k,  k  +  k^  X 


xt 


aoc. 


Hence,  at  time  t,  to  the  first  order 


Xj.  =  +  0^)  (y  +  Ky)  +  (S  +  6^)  +  k^(X+8^) 

....  •(  1  ) 


-  6  - 

CONFIDENTIAL 


C0NPID3NTLAL 


Technical  Ncte  Nc.  G-.W.566 

Now  if  and  are  the  true  position  and  velocity  at  a  time  T,  the  position 
at  any  other  time  T+t'  can  be  written 


T+t' 

+  V  ^ 

T 


(2) 


aiai  similarly 


T+t' 


r,T+t' 


=  X 


_  +  X  ^t'  +  f  5c 
•,T  r,T  J  1 


dt 


(3) 


B\it,  because  of  drifts  and  errors 


and 


•,T  “  ^  *  ®xT 


^rr  =  ^^^xT 


where  s^,  are  error  terms,  and  equation  (3)  becomes 


^r,T+t>  =  ^  ®xT  V‘  W’ 

T+t' 


f  [x+g^+(95^^+e^y)(y+gy)+(-0y^+0^^)(z+g^)+k^(3c+g^)]  dt^  .  (4) 

T 


Subtracting  (2)  from  (4) 


T+t' 


T+t' 


:,T+t'  ^+t'  “  ^xT  ^  ^xT*  *  J  ^  j  ^^zt^®xy^  (y+^) 


dt 


T 


T 


T+t'  T+t' 

+  I  (•‘9^yt"-®xz^  63,)  at' 


(5) 


PLCT/riting  this 


^xT"^  "^xT 


T+t'  T+t'  T+t' 

I't*  +  f  f  f 

T  T  T*^ 


=  X 


T+t*  T+t*  T+t* 

•,T+t'  •*  Vt'  ■  /  ®x  ■  /  ®xz^"^6z) 


(6) 
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and  putting  x  =  x  +  e,  to  the  first  order  this  is  equivalent  to 
S 


T+t' 


T+t' 


T+t' 


xT 


T  T 


~  ^r,T+t'  "  ^ 


T+t'  T+t'  T+t' 

g,T.t.  -  f  -  /  V/*'  -  /  * 

T  T  T 


(7) 


The  following  is  Chester's  approach.  Consider  a  sampling  period 
T-nv  to  T+nx  containing  (2n+l)  samples  at  eqxial  intervals  x.  For  each  sample 
there  is  an  equation  (?)•  Define  s  v  ,  "k.  ,  <f>  ,  <f>  as  those  constants  which 


X  A  jt  z  V  ^  2 

when  replacing  s^,  v^,  k^,  in  equations  (?)  reduce  \  e^,  over  all 

equations  to  a  minimum.  Then  if  e^,  is  random,  by  the  least  squares  criterion, 

these  will  be  the  'best'  estimates  of  the  values  of  k  6  -  and  6  For  a 

xT  ^zT  ^yT 

missile  constrained  to  fly  along  a  straight  line  near  the  vertical  this  method 
breaks  down.  Rewriting  equation  (6)  with  the  constants  we  have 


s  +  V 
X  X 


but  now 


T+t'  T+t'  T+t' 

j  ^2  J  (y-^Sy')^^^-<^y  j  =  f(t')  +  e^, 


T 
T+t' 


(z  +  sj 


T+t' 


* 


remains  almost  constant  over  short  sampling  periods  and  similarly 


T+t' 


(y  +  8y) 


T+t' 


(x  +  g^)  dt' 


and  the  equation  approximates  to  the  form 

T+t' 


3  +  V  t'  +  (k  +  c. ^  -  o„<p  )  /  (x  +  g  )  dt^  =  f(t')  +  e., 

X  X  ^  X  1^z  2^y'  J  ^  ^x  ^  '  t' 
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where  and  are  constants.  Hence  for  motion  in  a  near  vertical  straight 

line  the  values  of  k  ,  ^  and  <f>  are  not  separable  and  the  acceleration  errors 

X  z  y 

can  only  be  determined  if  the  drift  is  known. 

Provision  has  been  made  in  the  Black  Knight  inertia  table  experiment  to 
measure  drift.  Hence  equation  (6)  may  be  expressed  in  the  form 


T+t‘ 


T+t' 


""  ^r,T+t'  ■  ^g 


.T+t'  ■  / 


8xt“ 


T+t* 


T+t‘ 


(6) 


where  the  quantities  on  the  right  hand  side  are  measurable  and  can  be  sampled 
at  equal  time  intervals.  Replacing  the  integrals  by  the  appropriate  summations 
we  have 

i  k 

lc=0  j=50 


^r,T+iv  "*  ^g,T+iT 


-  T 


i  k  i  k 

2  ^  V‘  2 


^®x,T+jT”  J  J  ^^^z,T+jT'*'^xy^^^T+jT'*'S'>T-<-jT^^ 
k=0  j=0  kxO  j=0 


i  k 


-  T 


^  ^y,T+0T^^xz^^^T+ JT  ®z,T+jv^  ^  ®iT 


(8) 


IteO  j=0 


Two  ways  are  open  from  here.  Either  k  .  can  be  replaced  by  a  k  ,  without 

physical  counterpart,  which  is  assumed  to  be  constant  throughout  any  one 
sampling  period  and  whose  valve  oompvrfced  over  all  samples  is  used  as  the  value 
of  k^j  at  the  centre  of  the  period,  or  it  may  be  split  into  its  components  k^^ 

k^^,  k2^  ...,  which  are  constant  throughout  flight,  and  the  trajectory  from 

launch  to  bum  out  considered  as  a  single  sampling  period.  Additional  informa¬ 
tion  is  then  available,  in  that  both  velocity  and  position  are  known  to  be  zero 
when  t'  =  -nT,  and  two  rigorous  equations  must  be  solved  with  the  modified 
normal  equations,  using  the  method  of  Lagrangian  undetermined  mviltipliers. 

,  are  determined  directly,  the  computing  is  very  much  simpler 

when  done  from  launch  as  s  and  v  are  then  zero  and  need  not  be  estimated.  But 


When  kQ^,  k,^ 
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if  the  sampling  rate  is  highithe  numbers  involved  are  immense,  and  it  is 


f  1 

possible  for  rounding  errors  in  the  calculation  of 


E'‘  -  ( E“’) 


to  invalidate  the  answers.  Using  the  first  method,  by  taking  a  succession 

of  sampling  periods  ttooughout  the  flight  and  fitting  a  polynomial 

ku  +  k.  +  ku  X  +  ,,,,  to  the  values  of  k  ^  for  the  centre  of 

ux  lx  acc.  Zx  acc.  xT  aoo. 

each  period,  the  components  of  accelerometer  error  may  be  estimated.  Where 

the  number  of  components  is  uncertain  this  form  has  the  advantage  that  the 

degree  of  the  polynomial  finally  fitted  is  easily  adjusted. 


5  ANALYSIS  0?  ERRORS  DUE  TO  ANCILLARY  INSTRUMENTATION  AND  MEfflOD  OF 

COMPUTING 


k.  ,  k_  in  terms  of  the  errors  in  k 
^x’  2x ^ 

By  definition 


5.1  Errors  in  k^. 


X  _  X  k  _  =  k^  +  k.x  „  +  k„x 

aco.jT  xT  Ox  ix  aco.,T  ^x  aco.,T 


for  each  T.  But  the  estimates  of  k^  found  in  the  experiment  will  not  be 
without  error.  Call  the  error  e^.  Then  if  there  are  n  sampling  periods,  there 
will  be  n  equations 


acc, ,1 


acc. ,2 


X 

_  acc,  ,n 


X 


acc. ,1 


X 


acc,, 2 


X 


acc . , n 


1x 


2x 


x1  1 


^x2  ®2 


k  +  £ 

xn  n 


where  ^xl  is  now  the  experimental  value, 
Writing  this  in  matrix  form  as 


A 


n 


1x 


2x 


=  K 


n 


then,  solving  by  least  squares  for  the  best  estimates  of  ^x*  ^x* 
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=  A*K 

n 


=  A*6K 
n 


follow,  and  the  variances  of  in  terms  of  the  variance  of  k^ 

can  be  calculated,  with  the  proviso  that  if  the  variance  of  k^  changes 

systematically  through  the  trajectory  it  will  be  necessary  to  weight  the 
equations. 

5«2  Sources  of  error  in  k 

_ X 

The  equations  (8)  for  chosen  values  of  T  may  be  written  in  matrix  form 
for  a  period  of  m(=2n+l)  samples  as 


The  error  equations 


A*A 


^1x 


A*A 


6k 

6k. 


Ox 

1x 


.«>=2X 


A,q,  =  b+  o+  d+  e 
mx3  3  m  m  m  m 


where  q^  is  the  vector 


xT 

'^xT 

^xT 


and 


m 


0  = 
m 


m 


m 


r,T+iT 


-J  m 


“^g,T+iT  “  "" 


i  k 


EE  -.-E. 


k=0  0=0 


i  k 


_i 


-“r  >  >  (*i^2,T+jT  *  ®xy^  ^^T+Jt  S'^T+Jt^ 

k=0  j=0 

_i  k 

^"^y^T+OT  ®XZ^  ^\+jT  ®Z,T+jT^j 

k=0  j=0 


The  normal  equations  are  then 


A*Aq_  =  A*[b  +  0  +  d  +  e] 

^3  *•  m  m  m  m-* 

and  the  most  probable  value  of  k^  is  the  third  element  of  q^  in  the  solution 
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q,  =  [A  ♦A]  ^A*  [b+o+d+e]  . 

^3  L  j  '•  m  m  m  m'' 


If  6q,  is  the  error  in  q,  due  to  small  errors  6b.6o,6d,6e  inb.c 
^3  3  m  m  m  m  m'  m 


d  ,  e 
m'  m 


6q,  =  [A*A]  ^  a*  [6b  +60  +  6d  +  6e  ] 

^3  '•  m  m  m  m-* 


since  the  elements  of  A  are  constant.  In  particular,  if  p*  is  the  last  row 
in  the  matrix  [A*A]  ^A* 


6k  „  =  p* 

xT  -^m 


6b  +60  +  6d  +  6e 

m  m  m  m 


Since  the  sources  of  the  errors  in  b  ,  0  ,  d  and  e  are  independent 

m  m'  m  m  ^ 


E(6k  =  E(p*  6b  +  E(p*  60  +  E(p* 

'  jS!  m  m  m  m 


6d  )^  +  E(p*  6e 
m  '■^m  m 


and  the  contributions  to  the  r.m.s.  error  in  k^  from  the  various  sources  may 
be  considered  separately. 

5*3  Errors  in  recording  the  data 

Details  of  errors  inherent  in  the  different  methods  of  recording  the 
data  are  given,  with  estimates  of  the  variances,  as  appendices.  In  the  main 
they  fall  into  3  classes 

(1)  constant, 

(2)  random, 

(3)  slowly  varying  cyclic. 

Consider  the  errors  from  so\iroe  ’b*  only.  If  the  vector  p*  has  elements 
Pi,  -n  <  i  <  n 

6k^  =  ^  PjL  (9) 
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E(6k;^)2  =  6b2  *2^^  6b^6bj) 


i  <  ;3 


1  j 


=  ^  E(6b2)  +  2  ^  ^  E(6b^6bp  . 


-i*—  j 


(1)  If  5b^  is  constant  and  equal  to  cr^ 


E(6k^) 


2  /  2 


+  2 


-i  ‘-'j 


=  (r.l 


(2)  If  Sb^  is  random  with  variance  cr^ 


E(6k^)‘ 


2  2 
Pi  ^b 


since 


E(6b^6bp  = 


0  . 


(3)  If  6b^  is  slowly  varying  and  cyclic,  let  6b^ 
^  in  repeated  trials  is  random. 


a  sin(«t^+^)  where 


Then 


E(6k^)^  =  ^  p^a^  E(sin^(«t^+9i))  +  2a^^  ^  p^p^(sln(cjt^+^)sin((iatj+^)) 


i“j 


=  ^  E(l-cos2(wt^+^i))+a^  ^  ^  p^p^E|cosa3(t^-t j)-cos[oo(t^+tj)+2^]^ 


13 


1  2  V  2 

=  2^  p^ 


+  a 


PiPj  ® 


-1  <  e  <  1 
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2 

Now  if  a  is  the  an^ilitude  of  the  oscillation  -g-a  is  the  variance.  So  for  a 
slow  cyclic  variation  with  random  phase 


E(6k^)' 


lies  between 


4-2 


Pi^j  j  "-b 


and 


4^2 


Pipj  j  "'b 


5. 3»'l  Errors  in  position  due  to  errors  in  recording  acceleration  and 
their  effect  on  k 

X 


Taking  a  general  case,  s\;ppose  that  at  time  t  an  acceleration  is  X  and 

T+t’ 

that  the  error  in  the  recorded  acceleration  is  e  then  j  exdt^  may  be 

V  .. 

written  as  /  e>:(t’-q)dq  and  putting  this  in  the  form  of  a  summation  for 


sampled  data  using  the  trapezoidal  rule 


T+t' 


••2  2 

ex  dt  ±  -u 


••  HI 

(t  ♦  \  j 


i=1 


where  -u  is  the  sampling  interval,  e^x^  =  >.  niT  =  t '  and  the  sign  is  the 

sign  of  m.  If  6k  is  the  error  in  k  due  to  errors  ex  only,  then  it  follows  from 
(9;  that  ^  ^ 


/Equation  (10) 
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•• 

J’-n  (“'o  T  *  '-1^-1  ♦  ••••  ♦  ®-(.rl)^-(n-l)) 

•• 

P.(^l)(M)e„  ■#+  (n-2)e.,  +  ....  +  s.(„.2)^.(„.2)) 


♦  P-S  (2'o  T  +  =-1  ==-1 


*P-1  ro2 


+  p  X  0 
•^o 


P1  t; 

•• 

•  • 

f.  ^ 


^  f’n  'y"=o  T  + 


nr-1  n-1 


S  X 

-n  -n 


D  X  0 

-n  J 
+  '^.(n-1)  ^-(n-1)  !P-n>‘  ■' 


+  c 


-(n-2)  -(nr2) 


+  e 


b  2 


+  e 


X  r- 
o 


b  2 


^-n  +  P-(n-1 )_ 

np^n  P-(nr*l)  "^  ****  ‘^^-2  ^-1 

■*’  Pn-1  +  •••*  +  ^2 


(10) 


e  „  X 
n-2 

n-2 

•  « 

e  .  X 

n-1 

'r>-1  I 

••  I 

e  X 

fp 

X  0  I 

n  n 

J 

be  Vifritten 


5k  =  t: 
e 


E,- 
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vrtiere 


q.-  =  X. 

i 


[(n+i)  +  [(n-l)+i]  P„(n-1)  +  ••••  1 


-  n  <  i  <0 


%  =  [f  (I’-n  *  ^  (p-(n-l)  +  P(n-1)  *  *  P+l'] 


(n-i)  +  (n-1-i)  +  .... 


0  <  i  <  n 


5# 4  Contributions  to  the  variance  of  k  from  different  sources 

X 

5*4*1  Bias  error  in  recorded  acceleration  used  to  calculate  x 


This  follows  from  the  previous  paragraph  if  x.  is  made  equal  to  unity 

••  2  ^ 

so  that  the  error  dt,  is  independent  of  x.  If  a*  is  the  expected  variance 

^  2  2 
of  the  bias  errors,  the  contribution  to  the  variance  of  5k  is  a*  \ 


or  between  or 


. . .  4  -2  )  Vo 

the  errctrs  are  random,  constant,  or  correlated. 


0 


according  to  v/hether 


5. 4*  2  Scale  errors  in  the  recorded  accelerometer  output 

This  is  the  case  in  para.  5. 3.  "I  when  x  is  replaced  by  >  and  e 

represents  the  scale  factor. 

5.4.3  Errors  in  position  from  ground  instrumentation 

These  are  independent  of  acceleration  and  can  be  expressed  directly  in 
terms  of  the  p  as  in  para.  5.3*  If  is  their  variance  the  contribution  is 


2  V'  2 

o-  Wi 
i 


where  the  errors  are  random 


.)■ 


for  constant  errors 


and  between 


correlated.  Errors  in  position  Tiiich  are  systematic  and  can  be  expressed  as 
a  linear  function  of  time  during  the  sanipling  period  are  immaterial.  For,  in 
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the  equations  (6)  fhom  v/hioh  is  calculated,  if  =  x:^^^,+Ax+Avt*, 


T+t' 


xT 


*  =  ^‘r.T.t'  ■  ='g,T+t'  * 


can  be  written 


T+t* 

(s^  +  Ax)  +  (v^  +  Av)t'  +  j  (x  +  g^)at' 

T 


=  '^r,T+t*-'^+t*  + 

showing  that  the  errors  Ax  and  Av  are  absorbed  by  s  „  and  v  „  in  the  least 
squares  solution, 

5.4*4  Errors  in  drift  angle  measurement 

Since  the  drift  angle  enters  the  equations  as  a  factor  multiplying 
acceleration,  the  method  of  ccarputing  the  error  contribution  to  E(6k 

follows  that  used  for  the  acceleration  scale  error.  But  v/hereas  for  all 
other  sources  it  has  been  assumed  that  there  is  an  observation  with  its 
associated  error  for  each  sarnpling  instant  the  method  of  measuring  drift 
restricts  the  data  to  one  sample  per  second.  If  a  faster  sampling  rate  is 
used  intermediate  drift  angles  must  be  found  by  interpolation  and  will  be 
correlated  in  groups.  A  siniple  example  to  show  the  method  of  con^^uting  is 


given  below,  where  n  =  4»  t  ^  and 


i  odd. 


Erom  equation  (II) 
6k 


^-4  (P..4  X  0) 
^-2)  .. 


2  X  P„4) 

+  Eg  x^2  (2p_^^  +  p^^) 


(^2  +  O 


+  - 2~~^  ^-1  ^^P-4  ^-5 

•  • 

+  %  ■^f'**-(P-4+P4)+3(p«3+P3)+2(p„n+P2)+(p^1+P^)J 
+  (^Pi^  +  2p^  +  P2) 

+  62  (^4  +  P3) 


(Eg  +  E  ) 
—  -2“  "  ' 
+  %  \(o) 


*3^4 
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This  is  reviritten  as 


6k. 


2  " 


+ 


+  e 


-2 


-x__- 

~2  ^-2,  ■*■  ^-2  ^^-4  ■*■  ^-3^ 


3^  " 


-r  (^P-4+^-3‘'^-2)]  +  "=0  _  ] 


+  •  •  •  • 


i*e* 


6k 


e 


i  < 


n 

2 


and  the  standard  deviation  of  62^^  ,  i»e*  of  the  observations,  can  still  be 

used,  but  with  different  weights. 

5. 4»  5  Accelerometer  misalignment  errors  ' 

Since  the  misalignment  v/ill  be  measured  before  launch  it  is  only  the 
uncertainty  in  the  measurements  v/hich  need  be  considered.  This  will  give 
constant  errors  on  any  one  trial  and  the  contribution  to  E(6k^)‘^  will  be  given 

by  (T^  "inhere  cr  is  the  r.m. s,  error  of  measurement  and  x  is  replaced 

by  y  for  0  and  by  z  for  0  . 

acc.  xy  acc.  xz 

5. 4. 6  The  error  in  using  a  constant  'k*  throughout  the  interval  as  an 

estimate  of  the  value  at  the  centre  of  the  Interval 

Vhen  the  acceleration  can  be  expressed  as  a  linear  function  of  time; 
let  (x+g^)  =  a  +  bt  where  t  is  the  time  from  the  centre  of  the  interval,  T. 

(it  is  necessary  that  the  length  of  the  interval  is  such  that  for  maximum  t, 
bt/a  <  1.)  Then 


and 


xT 


Ox  ,  , 

- +  k.  +  k„  a 

a  1x  2x 


at  the  centre  of  the  interval.  In  each  equation  (6)  the  true  value  of  k 
can  be  written  as  ^,T+t 


^Qx  A  bt  J2.  >11  1  V4. 

—  “  T  "  (aj  "  “  ••••  )  +  N  +  V  ^2^^ 
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Thus  the  error  introduced  into  each  equation  by  taking  k  „  .  at  the 
central  value  is  ’ 


I  ^  (a+bt)  +  k2^bt(a+bt)J  dt^ 


T+t' 


=  k, 


[  -- 
Ox  /  a 

T** 


T+t* 

t  dt^  +  kn„  f  bt(a+bt)  dt^  . 
T 


^2x 


Since  kp^,  although  unknoi,7n,  are  oonstants  the  variance  of  the  error  in 

k^  is  given  by 


.-1 


vAiere  the  p^’s  are  the  elements  of  the  last  row  of  the  matrix  (A* A)  A*N 

T+t' 


•where  N  is  the  vector 


of  (A*A)*’^  A*M  ■where  M  is  the  veotor 


T' 


-  -  t  dt^ 

a 


and  the  q^'s  the  corresponding  terms 


..  T+t' 

Jbt(a+bt)dt^^  . 


6  ERROR  OdiPUTATIQN  FOR  AN  AEPROERIATE  BLACK  KNIGHT  TRAJECTCRY 


All  the  relationships  vAich  have  been  developed  for  the  X  co-ordinate 
apply  equally  well  to  Y  and  Z.  The  error  equations  have  been  programmed  for 
Pegasus  and  used  to  estimate  the  accuracy  vriL'th  ■which  it  might  be  hoped  to 
assess  "the  oomponent  errors  in  -the  accelerometers  in  the  inertia  table  ■txrials, 

A  •trajectory,  assumed  to  be  a  s^trai^t  line  at  2°  to  the  vertical  v/as  ocoputed 
against  time  from  launch  from  typical  B.K»  accelerations.  Approximations 
to  the  acoeleration  fiinotion  were  also  made  by  fitting  s^trai^t  lines  oen^tred 

on  10  secs,  30  secs . to  I30  secs  and  the  simpler  form  \ised  ■wi^thin  any 

one  sampling  period  for  the  method  of  computation  \7here  k  „  was  assumed 
constant. 

The  tables  1a  and  1b  ■which  follov/  show  the  ■weights  to  be  applied  to  the 
variances  of  errors  from  trials  instrumentation  for  different  sampling  rates 
and  sampling  periods  and  the  corresponding  E(6k)^.  Tables  2a  and  2b  show  the 
variances  of  the  errors  in  k^,  k^ ,  kg  in  terms  of  that  of  k  and  also  those  of 

k„  and  k.  viien  k_  is  absent. 

0  1  2 

The  variances  of  k^  and  k^ ,  and  k^,  k^  and  kg  calculated  directly  over 

the  -vThole  flight  are  given  in  tables  3a  and  3t>.  These  "were  oocputed  by 
integrating  from  launch  instead  of  from  the  centre  cf  the  "thrusting  period. 

^en  the  data  are  considered  as  a  whole  the  question  of  weighting  eqxiations  (8) 
beccanes  more  important.  But  the  choice  of  weights  is  difficult  as  ■the  differences 
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between  the  positions  from  ground  instrvimentation  and  those  estimated  by 
integrating  the  accelerometer  outputs  are  due  to  the  ccnibination  of  various 
errors  requiring  different  weighting  factors.  For  the  Black  Knight  experiment 
the  predominant  cause  of  error  in  kg^  was  found  from  the  earlier 

calculations  to  lie  in  the  ground  instrumentation.  For  this  only,  extra 
values  are  given  in  table  for  the  z  co-ordinate  showing  the  reduction  in 
the  variances  of  the  k*s  which  would  resxilt  from  weighting  the  equations, 
were  this  oonsidered  feasible.  The  question  also  arises,  7i/hat  errors  would 
be  made  in  k,^  and  kj  if  were  not  zero,  but  ignored?  These  are  shown  in 
tables  4.  1  ^ 


6. 1  Results  for  the  Z  accelerometer 

Calculation  of  k^,  followed  by  k^^,  k^  kg^.  By  taking  a  period  of 

30  seconds  and  sampling  at  second  intervals  the  r.m.  s.  error  in  k^  oan 

be  reduced  to  the  order  of  2  x  10  ^  and  with  10  such  periods  and  no  k^  term 

-2  -if 

present,  the  error  in  kQ^  is  of  the  order  of  1«3  x  10  g  and  in  k^^  1»3  x  10  . 

It  is  unlikely  that  these  values  could  be  inproved  much  by  increasing  either 
the  length  or  number  of  sampling  periods  as  the  xmture  of  the  instrumentation 
errors  is  such  that  correlation  would  arise  between  intermingled  samples.  So 
without  kg^  the  accuraoy  of  determination  of  this  method  is  hardly 

good  enough  and  for  k^^  it  is  not  good  enough.  The  plot  of  k^  z^  versus  z^ 

should  hov7ever  indicate  whether  or  not  a  kg^  term  is  to  be  expected  and  to 

this  extent  the  computation  wovild  be  of  seme  value.  With  kg^  present  the 

accuracy  is  too  lev/  to  assess  either  kQ^  or  k^  ^  satisfactorily. 


The  results  of  the  second  method,  whereby  k^^,  k^  ^  and  kg^  are  calculated 

directly  by  considering  the  whole  time  to  burn  out  as  a  single  sampling  period 

are  shown  in  tables  3a  sampling  intervals  of  5  seconds,  1  second, 

and  0.2  second.  For  the  z  accelerometer  this  appears  to  be  the  better  method 

but  the  accuraoy  is  still  too  low  to  meet  the  reqvdrements  for  either  k_ 

Uz 

or  k.  when  k„  is  determined. 

I  z 

6,2  Re  stilts  for  the  Y  accelerometer 


Using  the  first  method,  the  errors  in  k  are,  at  best  of  the  order  of 

-2  ^  — S 

1x10  and  the  method  fails  for  k.  and  k„  .  But  the  error  in  k^  ,  3  x  10  "^g 

1y  2y  Oy*  “ 

for  ten  thirty  second  sampling  periods,  suggests  that  it  should  be  possible  to 

detect  a  zero  error  in  the  Y  accelerometer  and  that  for  this  axis  the  first 

method  is  better.  For  Y,  using  the  second  method  the  error  introduced  in 

k^y  if  k^y  is  Ignored  is  of  the  order  of  0,7  k^y*  Hence  if  the  results  show 

%  to  be  small  the  errors  in  k^^.  in  table  3b  might  be  reduced  by  ignoring  k^^ 

as  well  as  kg^,  and  solving  for  k^^  only. 


7  DISCUSSION  OF  RESULTS 


In  all  the  computations  made  for  z  the  ground  instrumentation  random 
errors  were  predominant.  These  were  finally  weighted  as  for  the  radar,  on 
the  assumption  that  the  F.P.S.  l6  would  have  random  errors  of  0,0I!^  mils  in 
angle  and  8,5  ft  in  range,  giving  errors  in  the  z  oo^'crdinate  varying  from  k  ft 
at  launch  to  12  at  burnout.  Those  from  ballistio  cameras  would  be  apprwcimately 
half  these  but,  taking  a  sampling  interval  of  5  secs,  the  maximum  flash  rate 
now  being  quoted,  the  results  in  tables  3a  show  that  the  cameras  alone  will 
not  provide  the  answers. 
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Not  much  is  knovm  about  the  errors  in  the  radar  and  the  valxies  quoted 
may  v/ell  be  optimistic.  In  the  method  of  solution  ■which  estimates  errors  in 
position  and  velocity  at  the  centre  of  a  sanpling  inter'val,  any  systematic 
errors  "which  can  be  expressed  as  a  linear  function  of  time  are  absorbed  into 
these  'tv70  quantities  and  only  those  of  higher  order,  which  are  likely  to  be 
very  small,  ■will  contribute  to  the  acceleration  error  conponents.  But  in 
the  second  method  isiiere  the  burning  period  is  oonsidered  as  a  whole  and  the 
integration  is  done  from  lavinch,  systematic  errors  in  grotxnd  ins^trumentation 
can  be  serious.  Their  contributions  as  shown  in  tables  3a  and  3'b  were 
finally  computed  using  only  the  residuals  after  the  ballistic  camera  obsejrva- 
tions  had  been  used  to  remove  as  far  as  possible  the  radar  biases. 

Errors  in  the  estimation  of  drift,  "which  are  negligible  for  "the  z 

acceleraneter  calculations  are  as  important  as  those  from  the  ground  ins'trvimen- 

tation  in  the  calculation  of  k^  ,  k.  ,  k_  ,  and  because  the  tele"vlsion  cameras 

Oy'  1y'  2y* 

only  record  once  a  second,  their  oon"tributions  cannot  be  reduced  by  teiking 
fixes  more  often  than  this  since  it  "will  be  necessary  to  use  interpolated 
values. 

8  CONCLUSION 


If  the  accelerometers  to  be  tested  in  the  Black  Kni^t  inertia  table 
■trials  are  good  to  the  accuracy  specified  i.e,  2x10“^g  bias  and  2x10"^  scale 
error,  with  a  possible  non-linearity  in  the  soale  error,  the  method  of 
assessing  accelerometer  accuraoy  by  ccmparing  'the  integrated  output  with 
position  from  ground  ins'trumerrtation  may  prove  to  be  unsat isfao"tory  "with  the 
facilities  at  present  available. 


ACKNCr:7LEDGEMENT 


Thanks  are  due  to  UEr.  R.H.  Gooding  for  his  valuable  critioism  and  advice, 
and  for  checking  the  ocmputing. 


LIST  OP  REPERENCES 


Ref.  No. 


Author 


Chester 


2 

3 


Title,  etc. 

Notes  on  fli^t  test  data  reduotion. 

Chester  Ins"trumentation  Laboratory, 

MIT  Report  R-142. 

Letter  from  Miss.  Y/hitehead,  W.R.E. ,  dated 
22nd  February,  1960,  Ref,545V^7/2 

Evalua-tLon  and  analysis  of  radar  performance 
Ins'brumentation  Radar  AN/IPS  16(XIV-1)  1957» 

[g.\7.  Ref,  8955A] 


ATTACHED 

Appendices  1~5»  Tables  1-5>  Brgs.  GW/5p/9954  to  9956  Detachable  Abs'tract  Cards 
ADVANCE  DIgPRIEDTION;- 


M.O.A. 


R.A.E. 


Chief  Scientist 

CGRL 

DG"3M 

dgu?d(ta) 
G¥(G<SC)5  80 

TIL  90 


Director 

dd(e) 

dd(a) 

Trials  2 


Aberporth  2 
Patents 
Library 
Radio 


Arm 

lAP  3 

Instn 

Maths 


( 


-  21  - 

CONPD'JilNTIAL 


cowinmsiAL 


Technical  Note  No.  G.W.  5^6 


APPENDIX  1 

ERRORS  BETY/EEN  ACCELERaiETER  OUTFDT  AND  TEEJMETERED  DIGITAL  ACCELERATION 


Source 

(1)  Accelerometer  voltage 

supply; 

(2)  Anplifier: 


(3)  A-D  converter 

Reference  and  bias 
voltages: 
Resistors; 

Canparator ; 


Errors 


tenperature  variation 
vibration 

drift 

non-linearity 
tenperature  'j 
fluctuation  |  negligible 
vibration  J 


vibration 

tenperature  variation 
tenperature  variation  (negligible) 
time  variation 
rounding  errors 

uncertainty  in  round  off  level 
time  shift  in  round  off  level 


Random 

Description 

cr 

Comparator 

round  off 

1,2x10**^  ft/sec^ 

Con5)arator 

uncertainty  in  drift  of 
round  off  level 

2.2p<10“^  ft/sec^ 

Highly  correlated 

Amplifier 

zero  error 

1x10  ^  ft/sec^ 

Bias  voltage 

tenperature  variation 
(1  cycle/min) 
vibration 

1.1x10  ^  ft/sec^ 

<1 « 6x1 0“3  ft/ sec*^ 

Resistors 

variation  between  cali¬ 

For  B.K.  accelerations 

bration  and  firing 

ft/sec^  t=0-90  secs 
Mt^xlO"'^  "  t=90-120  secs 

15.2x10"^  «  t  >  120  secs 

Scale  errors 

Anplifier 

non-linearity 

6. 5x1 0"° 

Reference 

voltage 

tenperature  variation 
(1  cycle/min) 
vibration 

1.7x10“^ 

<  1x1 0“5 

Accelerometer 

voltage 

tenperature  variation 
vibration 

1,7x10“^ 

<  1x1 0“^ 

Constant 

Con^jarator 

uncertainty  in  roxmd  off 
level,  <  1/10  quantum 

2»4x10*"^  ft/sec^ 
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Bias  errors 


cr 


Random 


2.7  X  10**^  ft/sec^ 


Highly  correlated 


B.K.j  4x10^^  "  t  =  90-120  secs 

'<.5.5  X  "  t  >  120  secs 


Constant 


2.4  X  10“^  ft/sec^ 


Scale  errors 


Highly  correlated 


2.5  X  10 


NOTES;  (l)  If  the  errors  in  accelerometer  voltage  and  an^lifier  are  omitted 


on  the  grounds  that  these  form  an  integral  part  of  the  accelerometer, 
the  difference  in  cr  for  bias  errors  is  negligible  and  for  scale  errors 
O' reduces  to  2  x  10"5, 

(2)  The  errors  from  different  soiarces  have  been  combined  on  the 
assunption  that  they  are  independent.  In  the  case  of  those  caused  by 
tenperature  variation  some  correlation  is  possible  and  this  could  cause 
a  small  variation  in  the  cr  for  scale  errors. 
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AEPEIODIX  2 

ERRORS  IN  DRIFT  MG-LES  IviEASUKED  BY  TELEVISION  CA1>IERAS 


1  ERRORS  IN  CMERA  MEASUREMENTS 


On  the  television  screen  the  lines  will  be 
of  arc  apart.  If  the  focal  length  of  the  camera 


approximately  two  minutes 
is  f  and 


z 


c 


f 


(see  Pig. 3) 


as  the  star  appears  to  move  across  the  field  of  vievir  of  the  camera  the  readings 
of  3  vri.ll  be  in  steps  and  successive  readings  Tri.ll  be  correlated,  Pigs. 4  and.  5. 

The  relation  between  the  errors  will  depend  on  the  drift  rate,  which  is 
unknown,  and  the  most  that  can  be  assumed  is  that  the  error  in  a  single  reading 
is  random  with  a  rectangular  distribution  and  that  errors  in  successive 
observations  are  correlated. 

Por  a  the  errors  may  be  assumed  random  and  normally  distributed  vriLth  no 
correlation  between  successive  values. 

1 . 1  Relation  bettveen  the  direction  cosines  of  the  line  of  sight  to  a  star 

referred  to  fixed  space  axes,  and  the  measurements  fran  a  television 

camera 


In  the  Inertial  table  experiment  the  camera  axes  Tri.ll  be  in  the  YZ  plane 
Inclined  at  45°  and  135°  to  01.  Fig, 6. 

Camera  1 


The  direction  cosines  of  the  camera  eixis  referred  to  0,XIZ  at  laiinch  are 
(O,  I/V2,  1/^2).  Let  (■6,m,n)  be  the  direction  cosines  of  the  line  of  sight  to 
a  star  referred  to  the  fixed  space  axes  and  (■&  ,ni  ,n  )  those  referred  to  the 
camera  axes. 

Then 


e  =  e 

c 


ra 


c 


m+n 

/2 


n 


-m+n 

^^2 


Hence 


■6  m+n  -in+n  ^  ^  „ 

—  =  -  =  -  =  constant  C 

^c  /2f  '/'2  z 

c 

-  24  - 

CONFIDENTIAL 


CONPID^^NTIiiL 


Technical  Note  No.  G.W.566 
iippendix  2 


■6  =  X  o 

c 

m+n  =  ■/2  f  c 

-m+n  a  V2  2  c 

c 

(f  +  Zq) 

n  =  - c 

V2 

f-2 

m  =  - 2.  Q 

V2 

But  •6?  +  =  1 


and 


The  field  of  view  of  the  camera  is  only  +7*5  and  the  movement  within  it  is 
expected  to  be  small.  Hence,  since  the  camera  errors  are  quoted  as  angular 
errors,  the  direction  cosines  may  justifiably  be  written 
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e.  = 


m. 


1  -  3. 


/2 


:Jl  +  a^+3^^ 


(1  -  |3^)(1  “  +  3^^  + 

/2 


n. 


1  +  3. 


1 


V2 


Jl  +a^+3^ 


(1  +  3^)(1  -  ^ 
V2 


where  the  suffix  denotes  camera  1 


and 


6^^  5a^ 


and 


=  a 


63, 


&n. 


6n. 


1 


/2 

63^ 

V2 


O' 

m 


n. 


In  the  same  way  it  may  be  shown  that  for  camera  2,  the  direction  cosines  of 
whose  axis  are  (o,  — 1/V2,  I/V2),  the  direction  cosines  of  a  star  are 


,1^  * 4* 

-  (1  -  3p) 

“2  =  I - -  ■ .  ■  ■ 

V2\/1  +  ag  +  32 


1 


+  OLr 
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and 


6^2 


1,2  Calculation  of  drift  angles  in  terms  of  ■£-^m^n^;  ^^2^2 


If  the  directions  of  the  camera  axes  remained  constant  relative  to 
space  €ixes  a  star  in  the  field  of  view  of  the  camera  would  appear  stationary. 
If  the  camera  axes  drift  the  star  will  appear  to  move.  Let 

the  con5)onents  of  drift  about  the  fixed  space  axes  OXIZ,  Then  since  ^  is 
small,  if  we  assume  the  camera  axes  to  remain  fixed  and  the  star  sight  line 
to  suffer  rotations  ""^yt*  "^zt*  apparent  movement  in  the  field  of 

view  of  the  camera  will  be  the  same.  If  •6’m*n'  are  the  direction  cosines  of 
the  new  star  sight  line  referred  to  the  fixed  ^aoe  axes,  -fein,  ■C’m'n'  will  be 
connected  by  the  relation 


F—  W 

1 

i 

m 

t 

m* 

1 

! 

k 

n’ 

which  may  be  written 


0 


-n 

•HU 


+n 


0 


-m 

= 

^yt 

m*ia 

0 

n’-n 

(12) 

(13) 

(14) 


2  2  2 

Because  of  the  condition  6  -m  +n  «  1  the  three  equations  (12),  (l3)>  (14)  are 
not  independent.  But  taking  (l2)  and  (13)  frcm  the  first  camera  and  the 
corresponding  ones  from  cambra  2,  there  are  four  equations  in  three  unknowns 
which  may  be  solved  by  least  squares  to  give  estimates  of  Since 

the  errors  in  do  not  have  equal  weights,  it  is  better  to  revert  to 

a  and  P  and  write  the  equations  in  the  form 
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(■1+3^)  “(i-3^) 


-(•1+P^)  0 


0  d+Pg)  “Og"'') 


!  “(l+Po)  0 


-/ 2a^ 


'^xt 


yt 


‘^zt 


H  2  .2 

1  +a^ 

- 2~  2  ‘  ■^2a^'  -V2a^ 


v1+a^ 


■l+a^  +  3? 

-L - L.  (1-p.)  ..  (1..3  ) 


2  fl2 

+ap  +  P 

- 5 - s-  (-V^a*)  +  -/^a^ 


/1+a^  +P2 


1  +a! 


1*4 


0^1)  “  Oo“l) 


(15) 


To  the  first  order 


2  2 

1  +  +  P![ 

2  5 

1  +  a*'^  +  PJ'^ 


=  1 


and  the  right  hand  side  becomes 


vdiich  can  be  -written 


0 

0 


V2a' 

1 

1 

t _ i 

-^1* 

* 

-/2<x^ 

+  ttg 

I 

I _ 

V2  0 
0  -1 


0 

0 


0 

0 

-/2 

0 


r«  - 


^1* 


1 


Pl 


rl  - 


P* 

2 


-  P. 
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and  the  equation  (15)  1)60011108 


1 

72 

0 

c 

0 

0 

(1+3-,) 

-(1-e^) 

- J 

if 

J 

“  “1 

0 

-1 

c 

0 

-(1+^1) 

0 

/2a^ 

^yt 

^1*  - 

1 

-  72  0 

0 

0 

0 

(1+33) 

-(V) 

«2  “  «2 

0 

0 

c 

I  1 

-d+ej) 

0 

i 

^2  "  ^2  I 

!>• 

— 

0 

Ti  - 

72 

'^Xt'' 

oj  - 

1+3^ 

0 

-72a. 

1 

^yt 

3^-3, 

0 

1 

”72 

(1+32) 

4(P2-1) 

“2  ■  “2 

-(I+Pg) 

0 

W2a2 

^2  ”  ^2 

- 

where  the  a*s  and  3's  have  the  sajne  variance.  Expressing  these  equations  in 
matrix  form  as 


A,  ,  0,  =  1), 
4^3  3  4 


and  solving  by  least  squares,  the  normal  equations  (emitting  the  matrix 
dimensions;  are  given  by 


and 


A*A  ^  =  An) 

0  =  (A*A)“''A*b 


Now  for  any  non  singular  sqxxare  matrix 

6(0  =  -0^600^ 

SO 

6(A*A)  =  -(a*  a)  (A*6A+6A*A)(A*A) 

and 

60  a  6(A’*'A)"”*A*b  +  (A*A)"’’6A*b  +  (a*A)"''  A* 5b 

=  -(A*A)“''A*6A(A*A5'Vb-(A*A)"^6A*A(A*A)"'Vb  +  (A*A)’‘‘’6A*b  +  (A*A)~‘’A*6b  . 
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But  As 


a  tinit  matrix  of  dimensions  4. 


Hence 


6^  =  -(4*4)“^  A*6A(A*A)”^A*b +  (A*A)'‘”'  A^Sb 

=  -  (A^A)-”'  A*bk  <f>  +  (A*A)“”' 

But  ^  is  knovm  to  be  small,  so  to  the  first  order 


b<f>  =  (A* a)"”'  a*  5b  , 


and  the  errors  in  <f>  depend  vpon  the  positions  of  the  stars  in  the  fields  of 
view  of  the  cameras,  as  well  as  on  the  reading  errors. 

Some  idea  of  the  variations  in  errors  in  drift  angles  due  to  star 
position  may  be  found  by  considering  images  at  the  centre  and  boundaries  of 
the  field  of  view. 

If  suffix  1  denotes  camera  1  and  2,  camera  2, variances  have  been 
calculated  for  stars  at  0^  with  Og,  and  A^  with  Ag,  Bg,  ^2*^2  ^2* 

(Fig. 7). 

2  2 

The  table  below  gives  the  weighting  factors  for  o  where  o'  is  the 
variance  of  the  camera  errors,  assimied  the  same  for  a  and  3, 


°1°2 

A^B2 

A^C2 

^®2 

^°2 

E(5^^)^ 

1.0 

0.78 

0,98 

1.05 

0.85 

0.90 

E(5^6y)^ 

2.0 

1.56 

1.97 

1.97 

1,56 

1.77 

2,0 

2.43 

1.94 

1.96 

2.65 

2.27 

E(5^^5^y) 

0 

0 

0 

0 

0 

0 

E(5^^69i^) 

0 

0 

0.04 

0.36 

0.43 

0.21 

0 

0 

0 

0 

0 

0 
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jgEENDIX  3 

ERRORS  IN  POSITION  PRCM  GEOUND  INSTRUMENTATION 


1  BALLISTIC  CAMERAS 

Three  Tflld  B,C«4  cameras  v/^ith  Aviotar  lenses  (f  =  210  mm)  lie  in  a 
ring  of  radius  of  the  order  of  30  n.m.  about  the  launching  area  and  four  with 
Aatrotar  lenses  (f  =310  mm)  in  an  outer  ring.  The  latter  are  much  ftirther 
away  eind  less  suitable  to  cover  the  B.K.  burning  period,  7ith  the  following 
ncmenclature 


A  P.5 

B  ‘Toomera  Tech.  Area 
C  Coondambo 

D  Maralinga 
E  Oodnadattor 
P  Bereford 
G  Ceduna 


and  assuming  an  r.nus.  error  of  one  seoond  of  arc  at  each  camera,  the 
following  errors  in  position  result  fran  weighted  least  squares  solutions 
in  three  dimensions  for  a  standard  Black  Knight  trajectory. 


] 


inner  ring 


outer  ring 


z 

13,000  ft 

53 

,000  ft 

339,000  ft 

Stations 

0* 

cr 

cr 

cr 

1 

cr 

cr 

cr 

0* 

cr 

X 

y 

z 

X 

y 

z 

X 

y 

z 

A,B,C 

0,7 

0.9 

0.5 

0.6 

0.9 

0.5 

1.2 

1.2 

2.5 

A,B 

0.9 

1.9 

0.6 

0.8 

1.4 

a6 

1.4 

1.5 

3.2 

D,E,P,G 

3.0 

4.3 

2.3 

1 3.0 

4.3 

2.3 

3.0 

4.0 

3.0 

Standard  deviations  in  feet  per  sec  of  arc. 

The  accuracy  generally  quoted  for  the  cameras  vdth  Aviator  lenses  is 
5  seconds  of  arc  random  error.  Very  little  experience  has  been  gained  by 
0.7.  Dept. ,  R.A.E.  in  the  use  of  the  cameras  on  trials.  One  camera  plate  which 
has  been  calibrated  here  shovired  an  r.m.  s.  error  over  all  star  points  of 
5  microns  5  seconds  of  arc)  when  over  100  stars  ivere  used  to  compute  the 
plate  constants.  If  hovrever  as  few  as  12  stars  \7ere  used  to  calibrate  it  was 
possible  to  introduce  a  bias  of  the  order  of  3  seconds  which  wxxild  vary  over 
the  plate. 

One  set  of  observations'^  has  so  far  been  received  from  W.R.E,  These 
were  angles  frcjn  the  three  inner  ring  cameras  used  for  the  B.K,  05  firing. 

The  results  showed  a  random  error  of  approximately  5  seconds  of  arc  and  in 
addition  a  bias  of  the  same  order,  which  was  not  constant. 

The  effect  of  errors  arising  in  the  svirvey  of  the  camera  sites  has  been 
investigated  and  found  to  be  small  compared  with  camera  errors. 
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The  san^jling  rate  of  the  hallistic  cameras  is  limited  hy  the  frequency 
of  the  electronic  flash  lanit.  This  is  quoted  as  one  every  5  seconds. 

2  P.P.S.  16  RADAR 


All  information  on  the  accuracy  of  the  P,P,S,16  radar  has  "been  taken 
from  the  R.C.A.  report  on  Evaluation  and  Analysis  of  Radar  Performance^. 

The  errors  appear  to  fall  into  three  groups, 

(1)  Large  errors  due  to  propagation  and  dynamic  lag  for  which  it 
shoTjild  he  possible  to  apply  corrections. 

(2)  Bias,  both  in  range  and  angle  ifiiiich  will  vary  from  trial  to  trial 
but  shoiald  remain  constant  in  £iny  one  firing.  The  effect  of  this  on  a 
rectanguleir  co-ordinate  is  approximately  linear  with  time  for  Black 
Knight  during  a  san5)ling  period  and  as  such  will  contribute  no  error  to 

k. 

(3)  '•Noise'*.  The  noise  is  not  truly  random  but  has  low  frequency 
carponents.  This  reduces  the  value  of  vising  a  high  sampling  rate. 

The  figures  deduced  for  bias  are  of  the  order  of  0,C7  mil  in  elevation 
and  8.6  ft  in  range  and  for  noise  O.Off  mils  in  single  and  8»5  ft  in  range. 

These  have  been  reached  by  a  combination  of  theory  and  the  results  of  calibra¬ 
tion  trials  on  relatively  short  range  targets,  with  the  longer  ranges  generally 
associated  with  low  elevations. 

The  maximum  saii5)ling  rate  is  60/sec  and  that  proposed  by  W.R.E. , 

40/sec. 
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AEPENDEC  4 

THE  USE  OF  BALLISTIC  CAJiERA  OBSERVAIICMS  TO 

REDUCE  THE  RADAR  BIAS  ERRORS 


1  ESTIMATION  OF  BIAS 

1,1  The  siting  of  the  B,C,4  cameras  and  the  E,P, S,  16  radar  at  TTocsnera  is 
shown  in  Pig.1.  During  the  thrxist  period  of  Black  Knight  best  coverage  is 
given  by  the  three  cameras  in  the  inner  ring  and  in  estimating  positional 
errors  from  optical  data  it  has  been  assumed  that  observations  are  available 
from  all  three  stations.  These  have  been  used  to  oonpute  radar  co-ordinates 
of  the  missile,  R,6,0  where  0  is  the  elevation  and  <j>  azimuth  measvcred  frcan 
the  y  eucis,  and  their  covariance  matrices  for  varloxis  positions  along  the 
trajectory,  with  the  assunption  that  the  cameras  have  a  random  angular  error 
of  5  seconds  and  no  bias. 

At  any  time  t  let  (R,6,0)  Fig,2  be  the  true  radar  co-ordinates, 

those  given  by  the  radar  and  (R  0  0  )  the  values  calculated  ft’cm  camera  data, 
^  ^  o  c  c 

Let  AR  be  the  systematic  error  in  radar  range  and  AR^^  the  random  error  and 
let  AR^^  be  the  corresponding  random  error  in  camera  range  at  time  t. 

Then 


R  4. 

rt 

+ 

11 

AR-  + 

“rt 

and 

4. 

ct 

=  R  + 

^ct 

Therefore 

AR 

=  R  4.  -  R  4. 

rt  ct 

♦  \ 

where  is  randan  if  we  assume  no  bias  in  the  camera  observations.  Taking  a 

ntmber  of  times  t  and  solving  the  resultant  equations  by  least  squares,  the 
best  value  for  AR  is  the  mean 

i.  e,  AR  =  - =■ 


(Spt  -«ct  ) 


and  similarly  for  A0  and  A0. 
1,2  Residual  bias  errors 


The  randan  errors  in  radar  range  are  virtually  independent  of  distance 
and  so  have  a  constant  variance  (X^, 


Hence 
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and 

°Ae  2 
n 


2 


The  errors  in 
errors  in  4R, 


AR,  AS  and  A^  are  now  correlated.  For  if  6(AR),  6(Ae) 
AS  due  to  random  errors  in  R^^,  6^^, 


are 


6(AR)  6(Ae)  =  -V  r  y  ®rt 
n 

Here  the  errors  in  radeir  are  independent  of  camera  errors  and  the  radar 
range  errors  are  assvaned  to  be  independent  of  radar  angular  errors  and  so 


Likewise 

(“.e  %)  =  jE^  (x  Vi) 

2  RESmiAL  SYSTEMATIC  ERRORS  IN  y  AND  z 


From  Pig.  2 


Ay  =  AR  cos  6  cos  ^  -  R  sin  0  cos  ^  Ae  -  R  cos  0  sin  <f)  h<t> 
Az  =  AR  sin  0  +  H  oos  0  A6  • 
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Although  the  systematic  errors  in  the  radar  vary  from  day  to  <3ay,  on 
any  one  trial  they  may  be  considered  to  be  constant.  Hence  if  are 

the  variances  of  the  errors  in  y  and  z  caused  by  residual  systematic  radar 
errors  after  correction  by  ballistic  cameras 


Ay 


22  2  222  2  222  2 

0"^  cos  0  cos  (f)  +  R  sin  6  cos  ^  R  cos  0  sin  <f) 


N  2  2 

**  ^  ^'^AR  °A0^  ^  ®  ®  ^  cos  0  cos  <P  sin  4> 


AR 


2 

R  sin  0  cos  0  sin  <f)  cos  ^ 


2  2  2  2  2  2 

=  O'atj  ®  +  O'aa  R  cos  0+  (c.^c.  .)  R  sin  20 

Az  AR  A0  AR  A0 


Using  the  covariance  matrices  for  (R^^  0^^  given  in  tables  5> 


■^AR 

2 

^A0 

2 

^L<f> 


°AR°A0 


°AR°A^  “ 


4.1  ft 


0. 9  X  10  ^  radians^ 


1,3  X  10**^*^  radians^ 


1.7  X  10  ^  ft  radians 
1.0  X  10  ^  ft  radians 


A0  A^ 


-10  ? 
*0.73  X  10  '  radians 


These  give  con5)uted  values  of  cr  .  increasing  from  1.5  at  t  =  35  secs  to 

2  o  2 

4.6  ft  at  t  =  135  secs  and  cr  ?  remaining  practically  constant  at  2,1  ft"^. 

yc 
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UPPEmU.  5 

ERRORS  CAUSED  BY  kg  TOEN  THERE  ARE  OSCILLATIONS  IN  THE  ACCELERATION  OF 
FREQUENCY  HIGH  COMPARED  WITH  THE  TETJWTPY  SAIMilNG  RATE 


Let 


X  =  f(t)  ^  °os(a)^t+0^) 


Error  due  to  kg  =  kg  x 


cos(a)^t+0^)^ 


=  k. 


(f(t))^  +  2f(t)  a^  oos(a)^t+^^)+  ^  ®'i  cos^(a)^t+0^) 


Y  ^  cos(a)^t+0^)cos(wjt+0j)J 


i  d 


Error  due  to  k.  =  k.x 
1  1 


=  k. 


cos(a)^t+0^)n 


Error  due  to  kQ  =  kQ  . 


Let  x^  be  the  recorded  acceleration,  since  this  is  the  mean  value  taken  over 


15 


sec, , 


/Equation 
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**15 

j  (x+e)dt 


J_ 

15 


Djhere  e  is  the  contribution  from  aoc.  errors, 


=  15 


1 

t+7=*  ^+Tk  1=+7c 

j  f(t)dt  j  ^  ^i  j 


t  i 


.  1 
*^15 


.  1 
*^15 


+k^  /  y^i  °os(a)^t+^^)dt  +  kg  j  [f (t) ]^dt+2k2  j  oos(cc|j^t+0^)  dt 


+k2  j  cos^(co^t+0^)dt+k2  j  ^  a^a^  cos(a)^t+?i^)cos(a)jt+?ipdt 


i  d 


Nov?  if  we  assime  that  cc^  is  large  con5)ared  with  the  telemetry  sainpling  rate 


t+ 


1 


J  cos(a)^t+0^)dt  j  cos(a)^t+0^)dt 


""15 


J_ 

15 


_1_ 

15 


£:  0 


and  since  are  randan 


El 


a^a^  cos(ccij,t+^^)  cos(a)jt+^j)  0 


i  d 
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Hence 

_  2 

=  fit)  +  kQ  +  f(t)+  kg  [f(t)]^  ^2  ^  • 


But  f(t)  changes  slowly,  so  In  1/15  sec 


-  f — 

[f (t)  r  - 


and  we  may  vTrite 


X 

r 


f(t)  +  kQ  +  k^  f(t)  +  kg 


Y?hen  there  is  no  error  In  the  accelerometer,  i«e,  k^  =  k^  =  kg  =  0, 
X  ts  f(t),  and  the  k’s  have  been  con^nited  on  the  assuiDption  that 


X 

r 


f(t)  +  k.  +  k,  f(t)  +  k. 


f(t)) 


) 


ignoring  the  term  kg 


If 


remains  constant  throu^out  the 


aj 


bvirning  period  kg  —  vdll  appear  as  an  error  in  kQ,  since  they  vTill  not 

be  separable.  If,  on  the  other  hand,  the  power  of  the  spectrum  changes  with 
time,  it  will  be  necessary  to  \ase  wei^ts  and  errors  vri.ll  appear  in  the 
estimates  of  k.|  and  kg  as  well  as  in  kQ. 


Estimation  of 


A  vibration  record  relevant  to  the  longitudinal  accelerometer  has  been 
analysed  using  SPADA  for  B.K,  03.  The  original  record  consisted  of  seven  two 
second  bvirsts,  giving  vibration  frequencies  of  20  o/s  vqpwards  but  the  analysis 
is  unreliable  below  100  c/s.  The  mean  square  value  of  the  output  remained 
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relatively  constant  except  for  one  burst  at  50  sec  after  launch  VThen  it 
increased  by  a  factor  of  three.  Ignoring  this,  and  taking  an  average  value 
of  spectral  density  from  the  curves  as  approximately  8  x  10"^gvc/s  and 
integrating  frcm  20400  cycles/seo  the  maximum  frequency  to  'vjhioh  the  Kearfott 


-1  2 

X  10  g  .  The  estimated  value  of  kg 
to  a  ccmbination  of  kg  and  oscillations 


3  X  10”^  X  2  X  lo”^  X  1*024  X  10^  6  X  10”^  ft/seo^ 


a. 


accelercmeters  are  sensitive,  — =—  3 

-7  ^ 

la  2  X  10  and  thus  the  error  in  k^  due 

in  acceleration  would  be  approximately 


and  this  is  negligible. 
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TABLE  2a 


Errors  in  k2^  in  terras  of 


Factors  raultiplying 

aZ 


No,  of  saii^jling 
periods 

5 

10 

20 

28 

k 

Oz 

53749 

18554 

7752 

5294 

Vf.27 

12.89 

4.90 

3.263 

k 

2z 

0.00157 

0.00033 

0.00011 

0.00007 

^Oz 

5692 

2547 

1240 

879 

^1z 

1.85 

0.7602 

1 

0.3488 

1 

0.2433 

TABLE  2b 


Errors  in  k  ,  k.  ,  k  in  terms  of  o' 

Oy  1y  2y  ky 


Factors  multiplying 


ky 


No.  of  sanjjling 
periods 

5 

10 

20 

28 

Qy 

0.36 

0.229 

0.130 

0.0978 

k 

1y 

2.53 

1.16 

0.558 

0.396 

k 

2y 

0.30 

O.08I 

0.030 

0.020 

k 

oy 

0.16 

0.10 

0.066 

0,45 

k 

1y 

0.57 

0.28 

0.14 

0.10 
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TABLE  h. 

Errors  due  to  ignoring 


Z 


k 

Oz 

k 

1  z 

(0t  i 

! 

•  ff  2 

I  ^ 

at 

X  =  5  secs 

X 

CO 

• 

I  k  X  10“''^ 

1 

3.1  X  10^ 

4  X  10"”'^ 

a  =  1  sec 

! 

3.3  X  10^ 

4  X  10-^'- 

2.8  X  10^ 

k  X  10"^'" 

approx, 
r.m.  s.  error 

1.2  X  10 

fps^ 

3.4  X  10“^ 

j 

Y 


^Oz 

k 

1  z 

(Ut 

(Ot 

CM  CM 

T  =  5  secs 

3.8  X  10° 

4  X  10“^ 

1.5  X  10^ 

4  X  10“^ 

-1  =  1  sec 

3.0  X  10° 

4  X  10“^ 

1.3  X  lo"* 

4x10“^ 

approx, 
r.nus.  error 

3.7  X  10“^  fps^ 

7.5  X  10"^ 
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TABLE  5 

Covariance  matrices  for  R.6,^  calculated  tislng  ballistic  cameras 

(Provided  by  Ifr,  R.H,  Gooding) 


J  o 

^c*^0c 

^0^00 

i 

^QG^(fC 

.5lc^^ic 

00  <po 

<J^ 

(po 

R  is  measTired  in  feet  and  0  ani  <p  in  mils. 


t  =  16.5  secs 

10,1 

-0.00436 

0. 0526 

-O.OOif36 

0. 000372 

-0.000002 

0.0526 

-0.000002 

0.00152  _ 

t  =  82  secs 

10.1 

-0. 01 05 

0, 0366  ; 

! 

-0. 01 05 

0. 000395 

-0.000026 

0. 0366 

-0.000026 

0.00123 

t  =  113»5  secs 

14.3 

0.0089 

0.01A4 

0.0089 

0.000438 

-0.000045 

0.0144- 

-0.000045 

0.00104  _ 

t  =  1A6  secs 

129 

0.132 

-0.0305 

0.132 

0.000439 

-0.000006 

-0.  0305 

-0.000006 

0. 00220  _ 
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FIG.  I. 
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FIG. 3  -7. 


FIG.  3 -7.  DIAGRAMS  SHOWING  NOTATION  USED 
IN  TELEVISION  CAMERA  CALCULATIONS. 
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